Theory of dissipationless Nernst effects. 
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We develop a theory of transverse thermoelectric (Peltier) conductivity, a^y, in finite magnetic field - this 
particular conductivity is often the most important contribution to the Nernst thermopower We demonstrate 
that axy of a free electron gas can be expressed purely and exactly as the entropy per carrier irrespective of 
temperature (which agrees with seminal Hall bar result of Girvin and Jonson). In two dimensions we prove the 
universality of this result in the presence of disorder which allows explicit demonstration of a number features 
of interest to experiments on graphene and other two-dimensional materials. We also exploit this relationship in 
the low field regime and to analyze the rich singularity structure in axy{B, T) in three dimensions; we discuss 
its possible experimental implications. 

PACS numbers; 



The magnetothermoelectric Nernst-Ettingshausen effect ^ 
has enjoyed renewed interest in recent years, first as a probe of 
superconducting fluctuations, and more generally, as a novel 
transport characterization of electronic correlations. Follow- 
ing the initial work on the cupratesfi'] strong magnetother- 
moelectricity was found in a variety of interesting materi- 
als. While precise theoretical treatment is lacking for most of 
these cases, phenomenological descriptions in terms of con- 
ventional weak-field quasiparticle transport theorylQ 0. El Si 
or effective classical hydrodynamic models ^ Isfl have been 
used with varied degree of success i^. 

In this letter we break from these earlier studies to treat 
finite field response directly, with no recourse to a low-field 
regime. Our chief accomplishment is the exact expression of 
the off-diagonal Peltier conductivity, a^ij, in terms of entropy 
of free fermionic carriers (see Eqs. l5ll6ll7land[T0]i. In two di- 
mensions we prove the universality of this expression (which 
also applies to Dirac fermions) in the presence of quenched 
disorder and compare it against available experimental data. 
In three dimensions we obtain, essentially in a closed form, 
the entire intricate singularity structure in a^y (as a function 
of magnetic field and temperature) inherited from the Lan- 
dau level spectrum which bears strong resemblance to ther- 
moelectric phenomenology of graphite[10]. We also exam- 
ine the weak field limit, B — > 0, where we predict a sim- 
ple, axy — —s/B, dependence on magnetic field and entropy 
density s. Quite generally, a^y is somewhat of a less stud- 
ied and, hence, poorly understood quantity, at least as com- 
pared to electrical conductivity or entropy. Thus, our basic re- 
sult directly linking a-^y and entropy (importantly, without in- 
voking the so called "entropy currents" used elsewhere in the 
literature HI) is useful both for simplifying computations 
but also for developing intuition. Even if only approximate in 
more realistic models (e.g. with inelastic processes ignored 
by us here), it gives some credence to empirical associations 
of strong Nernst signatures with singular rearrangements of 
electronic structure, e.g. phase transitions. 




FIG. 1: (color online) Hall "brick" in a confining potential V{x) 
(above) and its Fermi surfaces showing three occupied Landau bands 
(depicted using different colors). The spectrum is discrete along kz 
axis and continuous along x. The chiral surface states occupy non- 
flat portions of the Fermi surface and flow along brick's sides. 



small thermal gradient is determined by 

J = (T-E-a-VT, 
J'^ = Ta-E- K-VT 



(1) 
(2) 



CuiTent flow in the presence of weak electric field and a 



where J, jQ,r,E, cr, a, K are respectively the charge and en- 
ergy currents; (slightly modulated in space) temperature and 
electrical field strength; electrical, Peltier and heat conduc- 
tivity tensors, respectively. Peltier conductivity is usually ex- 
tracted from electrical conductivity and thermopower, S = 
(T^^ • a, measured in a zero-current configuration. 

For convenience we consider a particular. Hall "brick", 
sample shape (see Fig. 1) of finite extent along x~ and 
z— axis, although our results will be independent of this as- 
sumption. Ignoring spin and possible valley quantum numbers 
for the time being the free electron Hamiltonian in a Landau 
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V{x), (3) 



are the elec- 



where niz, m, V{x), B, — e and Ib = w 
tron masses in the z direction and x — y plane, confining 
potential along the x-axis, magnetic field along the z-axis, 
electron charge and magnetic length, respectively. The con- 
fining potential along the z-axis (not shown explicitly) is as- 
sumed to be featureless (hard walls), whereas one along the 
X-axis, V{x), has a more gradual rise as would be the case 
in graphite nano-ribbons 1121 li3|], for example. Consequently, 
the wave functions are standing waves along the z— axis, la- 
beled with a discrete set of fc^'s. If V{x) = the spec- 
trum of this problem consists of Landau levels dispersing with 

kz, en,fc, = iT^n + where w„ = |f (n + vq), with 

1^0 = 1/2 (but can be different generally|!l4]). Provided 
varies smoothly on scales set by the magnetic length, the spec- 
trum of the Hall brick can be reconstructed by an adiabatic 
shift of Landau bands at each kz and n, e.g. by absorbing 
V{x) into a spatially varying va- Thus, the Fermi sea of the 
Hall brick is a locus of points in the kz — ky plane where 
En.fc, < Ai, the chemical potential (see Fig. 1). It is bounded 
by a set of closed Fermi surfaces each made up of two flat 

segments in the bulk (with kp n = ^2rn'^ 
chiral surface states propagating in ±y directions on opposite 
sides of the brick . While detailed structure of these surface 
states will depend on surface properties, we generally expect 
them to be robustly present (due to the necessary chirality) 
and posess interesting transport anisotropics between direc- 
tions along the equilibrium current direction (±y) and trans- 
verse to it, by analogy^ with chiral metals in layered quan- 
tum Hall svstems il5Ul6l1 . The analogy is not precise, since 
our chiral states are two dimensional, completely delocalized 
along the surface and smoothly connected to the gapless bulk 
(layered Hall states typically have mobility gaps in the bulk 
and strongly one dimensional surface states). Their wave- 
functions are peaked nearest to the surface for the smallest 
kz ~ tt/Lz component of the Fermi momentum (see fig. 1), 
which may allow for their characterization via coherent sur- 
face tunnelling. Detailed discussion of these chiral surface 
states will be presented elsewhere 1 17|1 . 

Bulk transport properties of the Hall brick can be com- 
puted easily by applying quantum Hall edge formalism for 
each transverse mode kz, e.g. following seminal works 
of Halperin[18] and Girvin and Jonson[19] on the Hall bar 
(which itself is a very thin Hall brick with only a single 
kz = 'T^/Lz state occupied). We apply both electric field and 
thermal gradient along the x which induces net currents along 
the y direction (see Eqs. 1 and 2). Recall|18, 3, that al- 
though microscopic current distributions depend sensitively 
on the details of the confining potential bulk transport coef- 
ficients computed by integrating over these distributions are 
independent of edge specifics, as they must be. The three non- 
zero Hall electric/Peltier/thermal conductivities can be written 
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xCo, axy = and K^y = -^^2, 



spectively, with the help of Cq — Cq {kz)/ Lz and 



dm 

de 



(4) 



Here and elsewhere, the Fermi function is denoted by 



/(e) = 
constant 
holes a 



1/(1 



), with ks the Boltzmann 



These coefficients are for electrons, for 
reverse sign. Explicit expressions 



xy and l^xy 



for Cq can be obtained by a further variable change 
de -> df and a definition /„ = /(ftw„ + 1^ - fi), 
with a familiar result for co{kz) — X^n /« 
somewhat less familiar expressions for ci(fcz) = 



fn log /„ + (1 - /„) log(l - /„)] and C2(fc^) = 
+ /„ log^l//„ - 1) - log'(l - /„) - 2Li2(l - /„ 



where Li2(2;) is the polylogarithm function. Thus, we find 
that, up to an overall prefactor, a^y is the entropy per particle 
added over Landau bands and averaged over transverse 
modes: 

axy = [•^" -^'^ + (1 - fn) log(l - /„)] . 



(5) 

Eq. |5]is the basic observation from which the rest of this paper 
follows (thermal conductivity is left for future work| 17]). 

We start by examining Eq. |5] in the semiclassical regime, 
with weak scattering and 5 — > - this regime may be realized 
in some semiconductors ||20[ I21I1 . Rewriting Eq. |3]in terms of 
entropy per volume (i.e. entropy density), s, we obtain 



Ctx 



s 



(6) 



akin to axy — where n is the density of particles. Ab- 

sent interactions, a^y is even under charge conjugation, while 
axy is odd. At high temperature, for a non-degenerate gas 
with s = kBfi, these results are completely consistent with the 
purely classical drift in crossed electric and magnetic fields, 
with entropic forcejl^l, F = —ks^T, and effective electric 
field is E ^ F/e. At low temperature, the original B — 
Fermi surface of the problem is more or less intact with semi- 
classical wavepackets of quasiparticles precessing around it. 
Thus we expect an additional de Haas-van Alphen type mag- 
netooscillations on top of the 1/B envelope as well as con- 
comitant Shubnikov-de Hass traces in conductivity. Also, 
Boltzmann kinetics can be usedlstl to interpolate axy{B) to 
very low fields where scattering dominates and axy ^ B. 

The two dimensional limit of Eq. |5]is obtained by omitting 
the sum over kz modes and setting /„ — f{fujOn — /i) 

, 00 

= ^ ^ [/n log /„ + (1 - /„) l0g(l - /„)] (7) 

which shows a sequence of thermally broadened nearly sym- 
metric peaks of height — log 2 as a function of B or [i, 
centered at quantum critical points separating integer Hall 



3 



plateaux (where axy shows activated behavior). This structure 
was discovered and explored previously, albeit in the context 
of thermpower, by Girvin and Jonson lfl9l I23I1 (although these 
authors apparently made no identification with the entropy). 
We now provide a proof that an exact quantum critical result 
in Eq. |7]is, in fact, universal in the presence of disorder. 

For simplicity we consider a setting with generic, short- 
ranged quenched disorder, strong magnetic field and no 
boundaries. We shall make use of two well established facts 
about this system. First, strong localization in the pres- 
ence of magnetic field results in a sequence of integer Hall 
transitions 1 24], with axy{T = 0, /i) T'^n^i^'' ^ 
changing in a discrete series of sharp steps typically located 
at ^ = hwn near the centers of original Landau levels (and 
9{x < 0) = and unity otherwise). Second, the so-called 
"generalized" Mott formula 



J — c 



de 



;Co(Ai + e,r = 0) (8) 



de keT 

is an exact integral expression shown previously lEsl Esll to 
relate axtj{T = 0) and axy{T). Substituting the universal 
behavior of a^y {T = 0) into the generalized Mott formula we 
find Eq. |7] with an important adjustment that /„ should be 
evaluated at the true critical points, /i = fta)„ rather than the 
Landau levels of the problem without disorder. 

Simple as it is, this proof is important on a number 
of counts. Firstly, it corrects the findings of Jonson and 
Girvin|23], who studied the effects of quenched disorder us- 
ing self-consistent Born approximation and concluded that 
their earlier, clean limit expressions lfl9ll were not universal. 
Secondly, universality in the presence of disorder implies that 
any experimentally observed deviations must be due to ex- 
perimental uncertainties, electron-electron or electron-phonon 
interactions, disequilibration effects or other physics omitted 
thus far Recent experiments on graphene|26] found a surpris- 
ingly robust confirmation of Eq. [7] for most plateaux transi- 
tions even in the strongly disordered samples, which also hints 
at the broader universality of Eq. |7] e.g. with respect to band 
structure details. In fact, we have confirmed this universality 
theoretically for Dirac fermionsfl?]. Thirdly, the fact that the 
peaks are centered about the true critical points in the pres- 
ence of disorder provides a very important additional insight 
highlighting the difference between the density of states (as 
measured through magnetization 1 27]) and the entropy - while 
the former accounts for both localized and extended states, the 
entropy of localized states is strictly zero, hence activated be- 
havior of axy is preserved. Finally, the standard Mott formula 



(9) 



transition and grossly understimates it away, in the activated 
regime, as T ^ 0. 

Next, we consider the three dimensional regime of large L^, 
so that a quasi continuum of momenta kz appears below the 
Fermi level, so we can write simply 



ek 



h ^ 



dkz 
2tt 



[/„log/„ + (l-/„)log(l-/„)] 

(10) 

One low temperature (Sommerfeld) expansion of this ex- 
pression can be obtained by linearizing electron dispersion 
about each of the flat Fermi surfaces with Fermi velocity. 



eks TT^ ■r^ ksT 



t^xy 



h 3 2'KhvFn ' 

n— 



(11) 



where n^ax = W], the index of the highest occupied Landau 
level, is the integer part of ly = {/.i/hujc — i^o)- Defining ther- 
mal deBroglie lengths for motion along the z— axis, Atz — 
h/y/2kBTniz and perpendiculai^ to it Xt — h/ y^2kBTni we 
can rewrite 



ekB TT^ £b 



1 



h 3 Xtz^t -n 
^ekBjB_ -J.) - ^(1/2, 1 + - v)] 



h XtzXt 



(12) 



in terms of ^, the Hurwitz zeta functions (i = \/^^). This ap- 
proximate expression for axy vanishes linearly with T as T — > 
0, it also diverges at a set of critical fields Bcn = e{n+uo) 
defined by hun — M via one sided root-singularities, due to 
one dimensional van-Hove singularities appearing every time 
a Landau band empties, with kpti ^ 0. Clearly, Sommer- 
feld expansion breaks down around these critical points and 
we need to do better 

To that end we rewrite Eq. (TO] with 6„ = /3(/i — hojn) and 
fz = f{z^ - 6„), 



x-xy 



ekB 
hXT7. 



(13) 



Fib) ^ 



F{b) 



with Aa = 1 is violated in this dissipationless regime, 
essentially due to singularities of zero-temperature conduc- 
tivity. For concreteness we can use (|2| and the finite T 
expression for axy in the absence of disorder to obtain 
AaifijT) fa ■§iP\fujJn — /i| in the activated regime and 
Aa = 121og2/7r^ ~ 0.84 precisely at the transition - ob- 
serve that Mott formula slightly overestimates axy near the tion of a^ 



/oo 
dzfz log/, + (1 - /,) log(l - /,)(14) 
-OC 

— > 00 

(15) 

L3 



e for b —00 
2 at 6 = 
2.4 for 6 = 6,„i„ ? 



in terms of the crossover function F (which is everywhere 
positive), whose relevant regimes are listed above and dis- 
played in Fig. I2] Most importantly, in quantum critical 
regimes he\B ~ Bcn\/'m ^ kBT we predict ~ VT varia- 
-2(e//i)(fcB/AT.). 
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FIG. 2: Using effective mass parameters m = 0.05me, niz = SrUe 
and carrier density n = 10~^* /cm'' we plot tiie expected behavior 
of a^y across tiie quantum limit showing both ~ T and quantum 
critical, ~ VT (in the inset), behaviors. We take T = 0.5, 1, 2, 4[/C] 
for the ~ T plot, and T = .25, .125, .0625, .03125[ii'] in the in- 
set (identifyable by the fact that lower temperature gives sharper 
peaks). We compute a^y in units of ao — — -^^T' where A = 
. The quantum critical behavior of the scaling func- 

.^2m^fcB(lKclvin) 

tion F is characterized by a temperature and field independent value 
at the crossing point (signifying the crossing of the Landau band's 
bottom) and a temperature independent value at the peak T away 
from the crossing point. 



Experiments: The broad range of materials where our the- 
ory may be tested directly includes conventional degenerate 
and non-degerate semiconduct ors^ ill , semiconductor het- 

, bulk semimetals lfioi 



28 



esp. ones with pronounced 
Our discussion here will 



erostructures and graphener2( 
[sO], quantum critical points 
changes in the Hall number 13 1 1 
be limited to some very recent experiments on graphene and 
graphite - the relative abundance of data and features make 
these ideal testbeds. 

Several groups have recently examined thermoelectricity in 
graphenel 2g 28 , 2M with data broadly confirming the exis- 
tence of peaks in thermopower near plateaux transitions. In 



particular, analysis of the data to extract a^y reveals II26II a 
reasonable agreement with Eq. [T] although the experimen- 
tal values appear somewhat smaller than the theoretical pre- 
diction in the entire doping range studied. It is unclear at 
present whether the observed discrepancy is within experi- 
mental uncertainty or if additional theory of inelastic pro- 
cesses (e.g. electron-electron or electron-phonon interactions) 
is necessary. As the non-interacting theory predicts peaks 
whose heights are independent of temperature, with widths 
(in field or gate voltage) narrowing as T ^ it would be of 
particular interest to measure the temperature dependence of 

Our theory for three dimensions appears to capture gross 
experimental features in graphite|10| rather well. The over- 
all scale of the observed signal is reproduced by the effective 
mass parameters for two dominant electron and hole bands 
previously used|32] in graphite: the estimate of a^y was ob- 



tained from thermopower data in Ref. (ld\ by taking diagonal 
resistivity p{B) w 2 • lO^^B 57 • m per Tesla and neglect- 
ing axxPxy contribution. Most importantly, the experiments 
clearly demonstrate the existence of asymmetric singularities 
synchronous with Landau level emptying transitions. Well 
away from these transitions Uxy usually vanishes linearly with 
temperature, as expected by the simple Sommerfeld expan- 
sion of Eq. [TT] While showing a clear breakdown of ~ T scal- 
ing near the singularities, the existing data lfioll is insufficient 
to confirm the ~ Vt law discussed above, unfortunately. 

In conclusion, we identified carrier entropy with axy ex- 
actly, irrespective of temperature. We demonstrated that 
strong localization in two dimensions enforces this relation- 
ship, computed axy in three dimensions and examined avail- 
able experimental data in graphene and graphite. While there 
is considerable degree of qualitative and quantitative agree- 
ment between our theory and the experiments, much work 
remains to appreciate the full reach and importance of our 
observation. Clearly, some theoretical treatment of inelas- 
tic processes is needed - while it appears [13] that weak 
electron-electron interactions leave our result intact, phonon 
effects (e.g. drag) on the Nernst thermopower are less 
understood [20]. Careful direct measurements (or extractions) 
of axy will be of great use as well. Finally, strongly correlated 
regimes, both at zero[33] and finite [!34l[35[l temperatures may 
also exhibit a version of Q;2:j,-to-entropy correspondence. 
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